The rise of spin-flip transitions in the anomalous Hall effect of FePt alloy 
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We carry out ab initio calculations which demonstrate the importance of spin-flip transitions for 
the intrinsic anomalous Hall conductivity of ordered FePt alloys. We show the such transitions get 
enhanced by large spin-orbit coupling of Pt atoms, becoming negligible when Pt is replaced by lighter 
isoelectronic Pd. We find that spin-flip transitions in FePt originate not only from conventional band 
anticrossings at the Fermi level, but also from transitions between well-separated pairs of bands with 
similar dispersions. We also predict a strong anisotropy in the anomalous Hall conductivity of FePt, 
which comes from spin- flip transitions entirely, and investigate the influence of disorder on it. 



The intrinsic anomalous Hall effect (AHE) |l| and 
spin Hall effect (SHE) in solids arise from the op- 
posite anomalous velocities experienced by spin-up and 
spin-down electrons as they move through the spin-orbit- 
coupled bands under an applied electric field. In para- 
magnets, where the bands are normally spin-degenerate, 
these counter-propagating transverse currents result in a 
time- reversal (T) conserving pure spin current. In ferro- 
magncts, where the bands are split by the exchange in- 
teraction, the same process generates a net T-odd charge 
current. 

The above picture is intuitively appealing, and often 
leads to correct conclusions. However, it leaves out the 
fact that in the presence of spin-orbit coupling the spin 
projection along the quantization axis is not a good quan- 
tum number. This is a particularly subtle point regard- 
ing the SHE, as the proper definition of the spin current 
becomes problematic when spin is not a conserved quan- 
tity 0] . More generally, processes which do not conserve 
spin are known to play a role in phenomena such as spin 
relaxation Q and magnetocrystalline anisotropy |||. It 
is however usually assumed that such spin-flip processes 
can be safely ignored when studying transport. This 
viewpoint is supported by recent calculations of the in- 
trinsic AHE [g| and extrinsic SHE Q. 

In this work we use first-principles calculations to 
study the impact of spin-flip transitions on the intrinsic 
anomalous Hall conductivity (AHC) of FePt orderered 
alloys 0. We find that their effect is considerable, as 
they account for about one fifth of the net AHC. More 
importantly, the calculations reveal a clear experimental 
signature of the spin-flip AHC: as the magnetization is 
rotated from the uniaxial direction to the basal plane, the 
spin-flip contribution changes sign, leading to a factor-of- 
two reduction in the net AHC, while the spin-conserving 
part is almost perfectly isotropic. 

We identify two distinct mechanisms for the spin-flip 
transitions. The first involves spin-orbit-induced anti- 
crossings at the intersections of up- and down-spin Fermi- 



surface sheets (the locus of intersection forms loops in 
/c-spacc, which we shall refer to as hot loops, in analogy 
with the hot spots that have been discussed in connection 
with spin-relaxation Q ) ■ The second mechanism involves 
spin-orbit driven transitions between bands with similar 
dispersion which are split in energy across the Fermi level. 
We shall refer to them as ladder transitions. Both fea- 
tures occur at very low frequencies, of the order of the 
spin-orbit coupling strength. 

Let us briefly review the formalism for calculating the 
intrinsic AHC from first-principles. For a ferromagnet 
with the orthorhombic crystal structure and magnetiza- 
tion M along the z ([001]) axis, the intrinsic anomalous 
Hall conductivity (AHC) a z = o~xv is given by the A:-space 
integral of the Berry curvature [l|, [9| : 

e2ft T f ,, V (V'rik | V x | -0mk) (V'rnk | Vy | ^ nk ) 
^ JBZ ^ (emk-Snk) 2 

(1) 

In this expression vb n ^ and tjj m k are respectively the oc- 
cupied (o) and empty (e) spinor Bloch eigenstates of the 
crystal, v x and v y are components of the velocity oper- 
ator v, and the integral is over the Brillouin zone (BZ). 
When the direction of M is changed from the z-axis to 
the x-axis ([100]), the a x = o~ yz component of the conduc- 
tivity tensor should be calculated instead, by replacing 
V x — > v y and v y — >• v z in Eq. (1). 

The calculations were done using the approach of 
Ref. whereby the linear-response expression (1) is 
rewritten in the basis of Wannier functions spanning the 
occupied and low-lying empty states. In this way the 
infinite sums over bands are replaced by sums over the 
small number of Wannier-interpolated bands. The Wan- 
nier functions were generated with WANNIER90 jTIJ using 
the same parameters as in Ref. by post-processing 
first-principles calculations done u sing the Jiilich DFT 
FLAPW code FLEUR [H] (see Ref. [IJ for details). The 
unit cell contained two atoms in the LIq structure, with 
stacking along the [001]-direction and lattice constants 
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Total 


a n 


a n 


Act 11 Act* 


FePt 


[001] 
[100] 


818.1 
409.5 


576.6 
585.1 


133.4 
-183.9 


-8.5 317.3 


FePd 


[001] 
[100] 


135.1 
275.9 


108.4 
196.9 


28.4 
62.0 


-88.5 -33.6 



TABLE I: Values of the AHC for [001] (cr z ) and [100] (a x 
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directions of magnetization M in FePt and FePd. Act 
defined as the difference between no-flip (flip) parts of a z and 
a x . All values are in S/cm. 



5.14 a.u. and c = 7.15 a.u. 14 



The spin-orbit term in the Hamiltonian has the form 
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where £ is the spin-orbit coupling strength and h is 
the magnetization direction, which is taken as the spin- 
quantization axis. The first term on the r.h.s. of Eq. (|2|), 
which we will denote as LS^, preserves the spin of a pure- 
spin state tpnk, while the second term, LS^, flips it. The 
spin-flip (<t*) and spin-conserving (a^) parts of the AHC 
a are calculated from Eq. (1) after selectively turning off 
the LS* or LS n parts of the SOI Hamiltonian ©. We 
find that to a very good approximation they are additive, 
i.e., a ~ cr" + <r*. 

The importance of spin-flip transitions for the AHC 
of FePt can be seen by analyzing its dependence on the 
magnetization direction (sec Table 1). If only the spin- 
conserving term in Eq. ^ is kept, the resulting AHC 
<jff changes by less than 2% from an average value of 
about 580 S/cm as M is tilted from the i-axis to the 
x-axis. When the the spin-flip term is also included in 
the calculation, the AHC becomes highly anisotropic, de- 
creasing by a factor of two from [001] to [100]. Keeping 
only the spin-flip part of the SOI reveals that it is indeed 
responsible for the large anisotropy, as the resulting AHC 
changes by more than 300 S/cm, from a positive value 
along [001] to a negative value along [100]. Such large 
AHC anisotropy can occur in uniaxial crystals, and was 
previously found in hep Co [l5[, however, as opposed to 
FePt, in hep Co the anisotropy is caused for the most 
part by spin-conserving transitions. 

Such a large spin-flip contribution in FePt is rather un- 
expected, in view of the fact that in a perturbative expan- 
sion in powers of £ only the spin-conserving part of the 
SOI Hamiltonian ^ contributes to Eq. (1), with spin-flip 
appearing only at second order [la ]. It should however 
be kept in mind that because Pt is a heavy atom, the 
SOI cannot be treated as a small perturbation in FePt. 
Moreover, the AHC is very sensitive to near-degeneracies 
across the Fermi level Q , and therefore the above analy- 
sis of Ref. [l6| , which is based on non-degenerate pertur- 
bation theory, may not apply. 

The AHC can be resolved in energy by defining a cu- 



FIG. 1: (color online) Spin- flip (%) and spin-conserving (ft) 
cumulative contribution to the AHC from the spectrum above 
energy uj, A(to). Inset presents the energy-energy density of 
contributions to the flip- AHC, S(ei, £2), for M along (a) [001] 
and (b) [100]-axes (in 10 5 a.u./eV 2 ). Open circles indicate the 
values of AHC from Table 1. 



mulative AHC A(lj), which accumulates all transitions 
in Eq. (1) for which £ m k — £ T ik > ^ [H|. In the limit 
lu — > all interband transitions are accounted for, and 
A(lu = 0) equals the full AHC. The spin-conserving and 
spin-flip contributions to the cumulative AHC are plot- 
ted in Fig. Q] in the range < u) < 12 cV, for both Mp 
and M||i. While A^ (lu) remains largely isotropic over 
the entire energy range and decays rather slowly with lu 
up to 4—5 eV in energy, the A^(lu) contribution picks up 
only for lu below 1 eV and immediately becomes strongly 
anisotropic with decreasin g e nergy, displaying a charac- 
teristic bifurcation shape [15} . Thus, the anisotropy in 
the AHC arises from spin-flip transitions in the 0.5 eV en- 
ergy window around Ep. 

To get a further insigt into which kind of transitions 
is responsible for A^-(u>), we calculate the density of con- 
tributions to the AHC given by Eq. (1), £(£1,62), from 
the states with energies £1 < Ep and £2 > Ep. The 
overall integral J J S(£i, s-2)de\ds2 provides the value of 
cr*, while assuming the constraint £2 — £1 > w, this inte- 
gral gives the value of A^(lu). Density S can be used in 
combination with the cumulative AHC to obtain more in- 
formation on the energy structure of the Berry curvature 
and the anomalous Hall conductivity. 

The calculated density £ for M||z and M||i is pre- 
sented in Fig. 1(a) and (b), respectively. In these plots 
we can clearly see the contributions from the band anti- 
crossings of |- and f-bands along the hot loops in the BZ. 
They are given by blue dots around the origin e\ = £2 = 
and provide a negative contribution to the AHC for both 
magnetization directions. While for M||i the hot loops 
contribution dominates, for M||z a competing positive 
contribution to the AHC can be clearly seen in Fig. 1(a). 
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Fe tot 


F# 


Fe n 


pt tot 




Pt tJ. 


[001] 
[100] 


-13.7 
210.0 


17.9 
253.6 


-26.8 
-37.5 


848.0 
65.0 


541.0 
425.7 


282.3 
-360.6 



TABLE II: Atomically resolved values of the AHC for [001] 
and [100] directions of magnetization M in FePt, decomposed 
into ft - an d ^-contributions. All values are in S/cm. 



It is given by scries of stripes ei — e\ ~ const, in the vicin- 
ity of the Fermi energy. By analyzing the band structure 
wc find that these transitions come from pairs of bands of 
different orbital character with similar dispersion around 
Ep (see inset in Fig. 2). Such ladder transitions, induced 
by SOI, provide a different source of the AHC as they do 
not require a band crossing at the Fermi energy, and oc- 
cur over large regions in energy and fc-space. In case 
of FePt with Mp their contribution is so large that it 
wins over the hot-loops part and determines the sign and 
magnitude of the flip- AHC. 

The spin-flip contribution to the AHC in FePt origi- 
nates from heavy Pt atoms. To demonstrate this we de- 
compose the spin-orbit part of the Hamiltonian in FePt 
in real space as 



H 



so 



Cf c L Fc • S 



£ Pt L pt -S, 



(3) 



where L M is the orbital angular momentum operator as- 
sociated with atom \i and £^ is the averaged over valence 
d-orbitals spin-orbit coupling strength. For an Fe atom 
in FePt the SOI strength £ Fc amounts to 0.06 eV, with 
corresponding value of £ pt = 0.54 cV. 

By using together representations of SOI according 
to (2) and (3), we perform a f?- and t-l-decomposition 
of the AHC coming separately from Fe and Pt atoms, 
presenting results in Table 2. It can be seen, that when 
we consider contribution to the AHC from only Fe atoms 
by setting £p t in (3) to zero, for Mp both f[- and 
AHC are very small, and while the Fe-driven cr^ remains 
also very small for the in-plane magnetization, dom- 
inates in this case. On the other hand, by analyzing the 
Pt-originated AHC (^Fe = 0), we observe that the spin- 
flip part is very large for both magnetization directions. 
For Mp the sum of Pt a 1 ^ and results in a large 
total AHC manifesting that for this magnetization direc- 
tion the AHE is driven by Pt atoms. For Mp, the Pt 
is of the same magnitude but of opposite sign to its jf- 
countcrpart, and both conductivities almost cancel — in 
this case the Hall current in FePt is mainly of Fe origin. 

When decreasing the SOI strength in FePt by substi- 
tuting Pt atoms with Pd atoms in FePd alloy, we see an 
essential decrease in the values of spin-flip AHC, as com- 
pared to the total AHC for both magnetization directions 
(Table 1), which is in correspondence to the perturbation 
theory arguments [lTj]- For FePd, the AHC anisotropy is 
mainly driven by Atr^, which is of the same sign and 
somewhat larger magnitude than Au^T j n FePt, while 




FIG. 2: (color online) Dependence of a~ (open circles) and 
a x (filled circles) in FePt alloy on the SOI strength inside the 
Pt atoms £pt with respect to its unsealed value £ Pt . Arrows 
indicate the spin-flip (%) or spin-conserving (ft) nature of the 
AHC anisotropy. Inset shows transitions in Eq. (1) leading to 
the ladder contribution to the flip-AHC. 



the total AHC anisotropy is opposite in sign to that in 
FePt, c.f. Table 1. In the following we make sure that 
such differences between the AHC of FePt and FePd are 
indeed caused by different SOI strength and not by other 
details of the electronic structure, such as Stoner param- 
eters, spread of d-functions etc. 

Without affecting SOI on Fe atoms, we scale down the 
£p t constant in FePt, Eq. (3), and perform self-consistent 
calculations of the a z and a x conductivities with the new 
corresponding SOI strength, £ Pt . The results of these cal- 
culations are presented in Fig. 2 as a function of £pt/£ Pt . 
In this figure we observe, that while the decay of a z is 
almost perfectly linear, the a x conductivity stays almost 
constant until the regime of £pt corresponding to the SOI 
strength of Pd atoms in FePd alloy £p d = 0.19 eV (indi- 
cated with shaded area in Fig. 2). This can be explained 
by the fact that for Mp the AHC comes from Fe atoms 
and thus it is not sensitive to £pt, however, the AHC 
for Mp is mostly Pt-originated (see Table 2), with £p t 
serving as effective SOI strength of the system in this 
case. 

At the value of £pt ~ £p t /2, a x starts dominating over 
o z , with values of AHC and its anisotropy qualitatively 
close to those in FePd, when £ Pt reaches £p d . We con- 
clude that the difference in value and sign of the AHE 
anisotropy between FePt and FePd alloys can indeed be 
attributed solely to the SOI strength of Pt and Pd atoms, 
and the crossover between [[- and ^-contributions with 
the SOI strength thus explains different sign of the AHC 
anisotropy in FePt and FePd, Fig. 2. Therefore, we sug- 
gest to use Pd-doped FePt in order to tune the effective 
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FIG. 3: (color online) Dependence of the cr z (open circles) 
and a x (filled circles) in FePt on quasiparticle damping rate 7, 
calculated within the Kubo-Streda formula. The area between 
ff- and "I4--AHC for two different magnetization directions is 
shaded in red and blue, respectively. 



SOI strength as well as spin-character and sign of the 
AHE anisotropy in these alloys. 

Finally, we investigate the influence of disorder on the 
AHC and its anisotropy in FePt. For this purpose, using 
the tight-binding formulation in terms of the Wannier 
functions, we employ the Kubo-Streda formula [HI for 
the AHC in the " constant 7 approximation" , assuming 
that the quasiparticle damping rate 7 is independent of 
the orbital [19| . The validity of this approximation for 
investigations of SHE and AHE in transition metals has 
been demonstrated 0, . The results of our calculations 
for a z (open circles) and a x (filled circles) as a function 
of disorder characterized by quasiparticle lifetime 1 /j are 
presented in Fig. 3, and allow for a simple explanation in 
terms of energy scales discussed previously. 

Upon increasing disorder, isotropic a" (red shaded 
area) stays practically constant until 7 of 0.1 eV, then 
decreases upon further increasing 7 and disappears at 7 
of 4—5 eV — a value, which characterizes the decay of the 
ff-cumulative AHC with energy, c.f. Fig. 2, and roughly 
corresponds to the characteristic band width in FePt. 
The t|-AHC (blue shaded area), on the other hand, de- 
cays much faster, and disappears at the value of 7 around 
0.5 eV corresponding to the width of in Fig. 2. This 
value can be traced back to the Pt SOI strength £p t (grey 
shaded area in Fig. 3), which emphasizes Pt origin of 
spin- flip contribution to the AHC. Upon 7 reaching £p t 
the interband coherence necessary for a build-up of t|- 
AHC is destroyed and cr* goes to zero. The relative 
robustness of with respect to 7, as compared to a^, 
underines the fact that the spin-flip transitions, living on 



a different energy scale, are much more sensitive to the 
degree of crystallinity, and arc affected stronger by disor- 
der. Overall, after adding up f[- and t|-AHC we observe 
that while the total a z (large open circles) monotonously 
decreases with 7 and drops significantly upon 7 reaching 
£ Pt , a x (larde filled circles) stays relatively constant in 
this range of disorder, with both conductivities vanishing 
at 7 of several eV. This qualitatively different behavior of 
<j z and a x upon increasing disorder brings us to a conclu- 
sion that up to a certain extent, the degree of disorder in 
FePt serves as the SOI strength £p t , c.f. Fig. 2, in accord 
to experimental findings (20j . 

To conclude, we predict a strong anisotropy of the in- 
trinsic AHC in FePt and FePd alloys. We show, that 
while in FePt the AHC anisotropy arises purely due to 
Pt-driven spin-flip transitions in a small energy window 
around Ep, upon descrcasing the SOI strength on Pt 
atoms, the sign of this anisotropy and its nature can be 
changed in Pd-containing alloys. We also demonstrate 
that in FePt the AHE comes from different types of atoms 
depending on the direction of the magnetization and that 
the degree of disorder in the samples of FePt can serve 
as an effective SOI strength of Pt atoms. 
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